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ABSTRACT 

We present spectroscopic and photometric data of the peculiar SN 2001gh, dis¬ 
covered by the Southern inTermediate Redshift ESO Supernova Search (STRESS) at a 
redshift z=0.16. SN 2001gh has relatively high luminosity at maximum (Mb = -18.55 
mag), while the light curve shows a broad peak. An early-time spectrum shows an 
almost featureless, blue continuum with a few weak and shallow P-Cygni lines that we 
attribute to Hei. Hei lines remain the only spectral features visible in a subsequent 
spectrum, obtained one month later. A remarkable property of SN 2001gh is the lack 
of significant spectral evolution over the temporal window of nearly one month sepa¬ 
rating the two spectra. In order to explain the properties of SN 2001gh, three powering 
mechanism are explored, including radioactive decays of a moderately large amount 
of ®®Ni, magnetar spin-down, and interaction of SN ejecta with circumstellar medium. 
We favour the latter scenario, with a SN Ib wrapped in a dense, circumstellar shell. 
The fact that no models provide an excellent fit with observations, confirms the trou¬ 
blesome interpretation of the nature of SN 200 Igh. A rate estimate for SN 2001gh-like 
event is also provided, confirming the intrinsic rarity of these objects. 

Key words: supernovae: general - supernovae: individual: SN 2001gh, SN 2002hy 


1 INTRODUCTION 

The most massive stars (larger than 25-30 Mq) are expected 
to lose their external hydrogen envelope in the final stages 
of their lives, becoming He-rich Wolf-Rayet stars (hereafter 
WRs) or, if they lose also their He mantle, C-O WRs. For 
decades, these stars were thought to explode as Type Ib 
or Ic supernovae (SNe), although there has been no direct 
detection of their progenitors to date. 

H depleted core-c ollapse (CC-) SNe display a variety 
of observed properties jTuratto. Benetti fc Pastorello|l2007l l 
depending on the stellar parameters (e.g. mass, radius, an¬ 
gular momentum) and the configuration of their circumstel¬ 
lar environment (chemical composition, density, geometrical 
distribution). Occasionally, SNe Ib/c display luminous light 


* E-mail: nancy.elias@oapd.inaf.it 


curve peaks and broad spectral features, which are indica¬ 
tive of high expansion velocity of the ejecta. They are la¬ 
belled as hypernovae, and became very popular since late 
90’ s for their connection with gamma-ray bursts (GRBs ; 


_GnJaina_etaL 19981 : IPatat et al.ll20oi] : Della Vallell2007l : 

iNomoto et al.ll2o'oW for reviews). There are also intermediate 
cases such as the Type Ib SN 2005bf, characterized by a pe¬ 
culiar, double-peaked light curve, that appear as transitional 
events between canon i cal SN Ib/c and gamma-ray bursters 
jArmpama et al]|2005l : iTominaea et al.llAosI : iFolatelli et al.l 

I2OO6II . 

Even more rare, a few extremely luminous SNe Ib/c 
were proposed to synthesize a large quantity of ^ ^Ni, e.g. the 
pecu liar Type Ic SN 1999as, for which 4 M© llPeng et al.l 
l200jl of ®®Ni were estimated. Because of their large in- 
ferred ^^Ni masses, objects like SN 1999 as and SN 2007bi 


dCal-Yam et al.ll2009l : lYoun^_et_^ 


be pair-instability SNe I Gal-Yam 2012h . although other 


[2010|) were proposed to 


© 2014 RAS 

































2 Elias-Rosa et al. 


plausible scenarios such as that driv en by the radioac¬ 
tive decay, were also suggested (e.g. iMoriva et al.l I2OI0I : 
iNicholl et al.|[2013ll . 

The operations of modern transient searches bring to 
the discovery of new types of objects that are challenging the 
classical SN explosion scen arios. Super-luminou s stripped 
envel ope/ H-poor SNe (e.g. Quimbv et ^ 1201 ll . iGal-YamI 


I 201 I and llnserra et al.ll2013l ) are an example of newly dis- 


s“^, slow-rising and bell-shaped light curves without a late¬ 
time tail declining following the ®®Co decay rate, and blue 
spectra. A few alternative scenarios have been proposed to 
explain these events, for example those where the SNe lumi¬ 
nosity has been powere d by the spin-down of a rapidly ro¬ 
tating young magnetar (iKasen fc Bildstenll2010ll . or by con¬ 
version of kinetic energy into radiation following the shock 
of shells with different velocities. The latter phenomena may 
result from many different evolutionary scenarios including 
non-disruptive ste llar explosions driven by pulsation al pair 
instability (see e.g. lWooslev. Blinnikov fc HegeI^2007^ . or by 
interaction of SN ejecta with a massive H- and He-poor cir- 
cumstellar medium (GSM, e.g. iBlinnikov fc So rokinal[2010l : 
IChevalier fc IrwinI I2OIII : iMoriva et al.l 2013l l or with more 
H-rich GSM, as in the case of the exceptional ly luminous 
CSS121015:004244-bl32827 (iBenetti et al.ll2014ll . The same 
mechanis m is invoked for the family of Type Ibn SNe 
(see e.g. iPastorello et al.l l2008l l. whose spectra are char¬ 
acterized by broad metal lines with superposed narrower 
He I features, and are believed to be core-collapse SNe ex¬ 
ploded within a He-rich GSM (e.g. iMatheson et al.l I2OOOI : 
IPastorello et al. 2008h. The prototype of this SN family is 
SN 2006jc ( Pastorello et al.ll2007l : iFolev et al.ll2007h . whose 
progenitor, likely a Wolf-Rayet star, was observed in out¬ 
burst shortly before the stellar core-collapse. 

Despite the rolling surveys have discovered thousands 
of transients and many new SN types, a mysterious object 
discovered more than a decade ago still remains an unique 
case. The peculiar SN 2001gh was discovered by the S'ortt/i- 
ern inTermediate Red shift ESO Supernova Search (STRESS, 
iBotticella et al.l I2OO8I : ESO stands for European Southern 
Observatory). The search was specifically designed to mea¬ 
sure the rate of Type la and CC SNe at intermediate red- 
shift. 

In the next section we will introduce the transient dis¬ 
covery. In section |3] we will present the spectroscopic (sec¬ 
tion EU and photometric data fsection |3.2II of SN 2001gh, 
while in section |4] we will briefly compare the case of 
SN 2001gh with other peculiar Type Ib SNe. In section [5] 
we will attempt to constrain via data modeling the explo¬ 
sion and ejecta parameters of SN 2001gh and, hence, char¬ 
acterize the progenitor star. Finally, in section [6l we will 
summarize our results. Additionally, an appendix has been 
included to illustrate previously unpublished data of a pecu¬ 
liar comparison object, the Type Ib SN 2002hy (Appendix 

Si. 


2 SN 2001GH 

SN 200Igh was discovered by lAltavilla et al.l (1200 ill on 
2001 November 12.19 UT (UT dates will be used through¬ 
out this paper). The object was in the STRESS field la¬ 



Figure 1. E-band image of SN 2001gh field taken with ESO 2.2 m 
-I- WFI on 2001 November 12. The stars used for the photometry 
calibration of the SN are labelled with numbers. 


belled J1888 JBotticella et al.ll2008ll at a = 00'‘57’"03f63, 
5 = -27°42'32'.'90 (J2000.0; see Figure [T}. SN 2001gh had 
a blue, featureless spectrum and was tenta tively classified 
a Ty pe II supernova at very early phases l|Altavilla et al.l 
2 OOII) , although this c lassification has later been revised 
Elias-Rosa et al.l ([200^. 

As we will see in the next section, the redshift estimated 
for this SN is z=0.16 ± 0.01. We will adopt a distance mod¬ 
ulus of 39.33 ± 0.10 mag, as derived from the Robertson- 
Walker metric considering z, Ho = 73 kms“^ Mpc“^, Qh = 
0.27 and Ha = 0.73. We also considere d a Galactic reddening 
correction of E{B — V) = 0.016 mag (ISchlaflv fc Finkbeined 
[2 oI 3). As there is no spectroscopic signature of interstellar 
Na ID at the host galaxy redshift, the contribution of the 
host galaxy to the total reddening is assumed to be negligi¬ 
ble. 

Unfortunately we cannot well constrain the explosion 
epoch of SN 2001gh. Given the strategy of STRESS, each 
field was observed on average every ~ 3 months. In partic¬ 
ular, our object was not detected on 2001 September 26, 
while the first detection is dated 2001 November 12, giv¬ 
ing us a window of approximately 47 days in the explosion 
time. Considering the epoch adopted for the B maximum 
light (see section |3^ and an ave rage value of 17.5 days for 
the rise of typical Type Ib SNe (IValentill2008l l. the explo¬ 
sion date for SN 2001gh could have happened just two days 
before its discovery, on JD = 2452225. 


3 OBSERVATIONAL ANALYSIS 
3.1 Spectroscopy 

Two spectra of SN 2001gh were obtained on 2001 November 
22.24 and 2001 December 20.16 using the ESO Very Large 
Telescope UT2, equipped with FORSl (and the grism/filter 
combination 300V-I-GG435; Figure EJ- Both observations 
were obtained at airmass ~ l'.'5 and with an average see¬ 
ing around or below I'.'O. The spectra were reduced follow- 
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ing standard iRAlj^ routines, including trimming, bias sub¬ 
traction, and flat-field corrections. The wavelength and flux 
calibrations of the one-dimensional spectra were performed 
using arc lamp and spectrophotometric standard stars spec¬ 
tra, respectively, obtained at the same night and with the 
same instrumental configuration as the SN spectra. A check 
on the flux calibration was performed using the available SN 
photometry and, for the second spectrum of SN 2001gh, a 
flux correction factor was applied to match the photometric 
data. 

The most notable feature in the first spectrum (ob¬ 
tained about 10 days after discovery) is the very blue contin¬ 
uum, with little evidence of spectral lines. A few weak and 
shallow P-Cygni profiles are observed at 3889, 4471, 5015 
and 5876 A (rest frame) that, in cont rast with the origi¬ 
nal identification (lAltavilla et al.l[2001^ . we now attribute 
to He I. The lines, in particular the A5876 feature, seems to 
have boxy profiles (see the insert panel in Figure ij. This 
is possibly an indication of the presence of an outer shell 
(see, for example, the boxy profile in SN 1993J; e.g. 
IPatat. Chugai fc Mazzalilll995f l. The second spectrum (ob¬ 
tained approximately one month later) shows a slightly red¬ 
der continuum, which is indicative of a temperature decline. 
Yet the Hel P-Cygni lines remain the only features visible 
in the spectrum, which does not show clear metal features 
commonly seen in ordinary CC-SNe. With the revised line 
identifications, and the lack of spectroscopic evolution be¬ 
tween the two epochs, we propose to reclassify this SN as 
a peculiar Type Ib. Considering the position of the emis¬ 
sion peaks of the He l lines, we constrain the SN redshift to 
z=0.16 ± 0.01. 


In Figure we show a comparison of these spectra with 
coeval spectra of the Type Ib SN 2008D (iMazzali et al.ll^OOSl : 
iModiaz et al.l 120141'). the super-lum inous stripped envelope 
SN 2010gx lPastorello_^tak 2010l l. and the normal SN II- 
P 1999em ([Leonard et al.ll2002l l. Neither the pronounced P- 
Cygni profile of Ha characterizing spectra of SNe H (6562.8 
A), nor the Fell around 5000 A of normal SNe Ib (SN 2008D) 
are unequivocally detected in the spectra of SN 2001gh. The 
Call H&K (3934, 3968 A) feature is instead visible in the 
spectra of the SN, though probably blended with Hel A3889. 
In contrast with the other SNe shown in Figure[21 SN 2001gh 
show very little spectral evolution, and the He l lines remain 
the only prominent features visible in the December 20th 
spectrum. 

The velocities, as derived from the position of the P- 
Cygni minimum of the most prominent Hel line, the A5876 
featur^, are ~ 6600 and 6000 kms ^ respectively, for the 
early and late spectra of SN 2001gh. While the line velocity 
determined after the maximum light is comparable with that 
of Type Ib SNe at coeval phases, the value estimated at the 
earlier epoch is relatively small, being the typical line veloc¬ 
ities of stripped-envelop e SNe at this phase ~ 10000 kms“^; 
e.g. ([Branch et al.ll2002 L This low velocity at peak may be 


® IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are oper¬ 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA), under cooperative agreement with the National Sci¬ 
ence Foundation (NSF). 

® The low S/N of our spectra does not allow us to properly mea¬ 
sure the velocity of other lines. 


an indication that SN 2001gh ejected more mass than typ¬ 
ical SNe Ib, and/or a lower kinetic energy is released in 
the explosion. Additionally, we have measured a black-body 
temperature by fitting the continuum of our spectra. It has 
decreased from Te// ~ 15000 K, as estimated in our first 
spectrum, to Te// 7400 K from the second one. These 
values are comparable with those of SN 2010gx or other 
super -luminous SNe at co eval epochs, as reported for exam¬ 
ple in [lnserra et al. lliH), but they are hi gh in comparison 
with th ose of normal Type Ib SNe (e.g. see [Stritzinger et al.l 
[ 2 OO 2 I or [Modiaz et aD[2009l) . 

3.2 Photometry 

The photometric images were reduced making use of the 
template subtraction technique to remove the possible con¬ 
tamination of the host galaxy light. Image subtr action was 
performed using the ISIS program l|Alar dl [200nl) . As refer¬ 
ence images (templates) of the host galaxy, we selected those 
taken about two years before SN 2001gh discovery on 1999 
December 12 and 1999 August 4 with the MPG/ESO 2.2 m 
telescope -|- Wide Field Imager (WFI) in V- and i?-bands, 
respectively. After geometrical and photometric registration 
of the two images (target image and template), and the 
degradation of the image with the best seeing to match the 
worst one, the template was subtracted from the target im¬ 
age. The instrumental magnitude of the SN was measured 
with the point-spread function (PSF) fitting technique us¬ 
ing the DAOPHOT package on the subtracted image. Refer¬ 
ence stars in the SN field were also measured using the iraf 
PSF fitting routine on the original image (see Table [T]). In 
order to calibrate the instrumental magnitudes to a stan¬ 
dard photometric system, we used the specific colour term 
equations for each instrumental configuration using average 
colour terms provided by the telescopes teams. 

Given the high redshift estimated for this SN (see Sect. 
13.111 . the observed V- and i?-bands are respectively close to 
the B- and F-bands in the rest frame. The K-correction for 
SN 2001gh was estimated by measuring the differences be¬ 
tween synthetic photometry in the observed and rest-frame 
SN spectra. Since we had only two spectra of SN 2001gh, 
we have assumed that the K-correction derived from the 
first spectrum applies to the early epochs, likewise, the K- 
correction derived from the second spectrum was applied for 
the photometry taken at the latest epoch. The K-correction 
for the intermediate epochs was estimated by interpolating 
the values derived from the two spectra. A K-correction of 
the order of a few tenths of a magnitude was estimated at 
all epochs. K-corrected magnitudes are used throughout this 
work (see Table [2j. Total uncertainties were computed tak¬ 
ing into account the photometric errors and the uncertainties 
in the K-correction. 

Although the photometric monitoring of SN 2001gh is 
rather sparse, it provides a good coverage of the SN evolu¬ 
tion during the photospheric phase (see Figure[3l top panel). 
Some upper limits were also measured before and after the 
B-band maximum, by introducing artificial point sources 
with progressively decreasing flux at the SN location. In or¬ 
der to estimate the epochs and magnitudes of the maximum 
light for the two bands, we fitted the two light curves with 
low-order polynomials, obtaining JD 2452242.5 ±1.3 and 
20.85 ± 0.10 mag for the B-band, and JD 2452247.1 ± 1.1 
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Figure 2. Comparison between optical spectra of SN 2001gh a few days before maximum light (a) and one month later (6) with those 
of SN Ib 2008D, the super-luminous stripped envelope SN 2010gx, and the SN II-P 1999em. The inset enlarges the region around the 
He I A5876 feature. Vertical solid lines mark the rest frame position of the main He I features, while the dotted, dashed, and dot-dashed 
lines give the rest wavelength positions of the main H and Fell lines, along with those of the Call H&K doublet. All spectra have been 
reddening and redshift corrected adopting the values derived in this paper and from the literature. 


and 20.70±0.20 mag for the V-band. We note that the broad 
light curve peaks are unusual, and not observed in normal 
Type Ib SNe. 

In Figure[3](bottom panel), the evolution of the intrinsic 
{B — V) colour of SN 2001gh in the rest frame is shown. The 
{B — V) colour of SN 2001gh is nearly zero at early epochs, 
and then becomes redder until it reaches (B—V) ~ 0.45 at 15 
days after the maximum light. At this epoch, the SN seems 
to start a blue-ward or a flat trend. Unfortunately our data 
stop, preventing us to check the true behaviour of this SN. 
As we can also see in the figure, the overall colour evolution 
of SN 2001gh is different from tho se of normal SNe Ib such 
as SN 2008D llModiaz et al.llioo^ l. Additionally, the colour 
of SN 2001gh is bluer than SN 2008D, which is consistent 
with the higher temperature shown in the spectra. 


Figure [4] shows the rest-frame U-band absolute light 
curve of SN 2001gh, compared with those of the same 
stripped envelope SN sa mple used for the spectral com- 
parison, viz. SNe 2008D llMazzali et al.l[2008l l. and 2010gx 
JPastorello et al.) 2010l'). and with that of the Type Ic 
SN 19941 I Richmond et al] Il996l l. If we assume for 
SN 2010gx a smooth photometric evolution, its light curve 
peak of SN 2001gh would be very broad, with a peak lumi¬ 
nosity that is significantly brighter than those of most CC- 


Table 1 . Magnitudes and associated errors of the stellar sequence 
used in the calibration process of SN 2001gh’s photometry. 


Star 

V 

(mag) 

R 

(mag) 

1 

17.11 ±0.01 

16.02 ±0.01 

2 

15.75 ±0.01 

15.02 ±0.02 

3 

18.73 ±0.01 

17.77 ±0.01 

4 

22.07 ±0.01 

20.86 ±0.01 

5 

19.81 ±0.02 

18.92 ±0.02 

6 

20.80 ±0.10 

19.45 ±0.01 


Note that the photometry of SN 2001gh was calibrated in V- 
and R-bands before being K corrected. 


SNe, although much fainter than that of the super-luminous 
SN 2010gx. 

Despite the incomplete coverage of the light curve, 
it seems that the SN luminosity declines quite slowly af¬ 
ter the maximum light. This behaviour could be indicative 
of additional powering mechanisms, such as recombination 
throughout a massive He envelope or even interaction with 
circumstellar material (see discussion in Section [5)). 
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Figure 4. F-band absolute light curve for SN 2001gh compared with those of SNe Ib 2008D, the super-luminous 2010gx (B-band 
absolute magnitude), and the Ic 19941. The magnitudes of SNe 2001gh have been corrected for the assumed reddening, E{B — V) of 
0.016 mag, and the adopted distance modulus of 39.33 ± 0.10 mag. Distances and extinction estimates for all other SNe have been 
adopted from the literature. Ages are relative to B maximum light in the rest frame. 


Table 2. K corrected B- and F-band photometry of SN 2001gh. 


UT Date 

JD 

2400000.0+ 

Phase^ 

(days) 

B 

(mag) 

V 

(mag) 

Instr.^ 

04/08/99 

51394.9 

-730.7 

- 

< 23.7 

WFI 

03/12/99 

51515.6 

-626.7 

< 24.5 

- 

WFI 

26/09/01 

52178.8 

-55.0 

< 23.0 

< 22.5 

WFI 

12/11/01 

52225.7 

-14.5 

21.66 + 0.10 

- 

WFI 

13/11/01 

52226.7 

-13.6 

- 

21.58 + 0.14 

WFI 

19/11/01 

52232.7 

-8.4 

21.07 + 0.10 

- 

WFI 

20/11/01 

52233.7 

-7.6 

- 

21.00 + 0.07 

WFI 

22/11/Olt 

52235.7 

-5.8 

20.93 + 0.07 

20.83 + 0.08 

FORSl 

09/12/01 

52252.6 

8.7 

21.06 + 0.10 

- 

WFI 

10/12/01 

52253.6 

9.6 

- 

20.70 + 0.31 

WFI 

16/12/01 

52259.7 

14.8 

21.25 + 0.10 

20.79 + 0.07 

EMMI 

20/12/Olt 

52263.7 

18.2 

21.23 + 0.07 

20.83 + 0.13 

FORSl 

06/06/02 

52431.9 

163.3 

- 

< 24.4 

FORS2 


^ Rest frame phases relative to the B maximum (JD = 
2452241.1 ± 0.2). 

2 WFI = ESO 2.2 m -h WFI, 0.238"pix-i; FORSl= ESO VLT 
UT2 -I- FORSl, 0.250"pix-i; EMMI = ESO NTT -I- EMMI, 
0.167"pix-i; FORS2= ESO VLT UTl -|- FORS2, 0.250"pix-i. 
f These dates are coincident with those of the spectroscopic 
observations. 


4 COMPARISON OF SN 2001GH WITH 
PECULIAR TYPE IB SNE. 


To place SN 2001gh in the context of the Type Ib SN zoo, 
we compare in Figure [5] our earliest spectrum of SN 2001gh 
with spectra of a sample of peculiar Ib SNe at similar 
epochs, viz. SNe 200 2bj (a peculiar, fast-evolving SN Ib; 
IPoznanski et al]l20ld) . 2002hy (a relatively narrow-lined SN 
Ib from the Padova-Asiago SN archive, cfr. appendix 0, 
2005bf (a Type Ib event with a double- peaked light curve; 
iFolatelli et al'l l2006l : lModiaz_et_^ ). a nd the unusual 
Type Ibn LSQ13ccw i Pastorello et al. 20151 ). Among them, 
an object that shares some spectroscopic similarity with 
SN 2001gh is SN 2002hy. The spectra of these two SNe 
are blue, with Hel lines showing evident P-Cygni profiles, 
such as A3889, A4471, A5015 and A5876. Weak Can H&K 
(AA3934,3968) features are also detected in the spectra of 
the two objects. 


The most remarkable difference is that the P-Cygni ab¬ 
sorption at around 5870 A is more blue-shifted in the spec¬ 
trum of SN 2002hy than in SN 2001gh. In fact, the velocity 
of the ejecta of SN 2002hy, as derived from the position of 
the P-Cygni minimum of the prominent Hel A5876 line, is 
8000 kms“^, higher than that estimated for SN 2001gh in 
the early epoch spectrum. In addition, SN 2002hy shows a 
cooler spectrum with a measured black-body temperature of 
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Figure 3. K-corrected B and V light curves of SN 2001gh (top) 
and comparison of the intrinsic (B — V) colour evolution of 
SN 2001gh and SN 2008D (bottom). The plotted V-band light 
curve has been shifted by 1 mag. Dotted vertical lines indicate 
the epoch when the SN 2001gh spectra were taken. Ages are rel¬ 
ative to B maximum light in the rest frame. 


^10000 K, versus ^15000 K estimated from the SN 2001gh 
one. On the other hand, the properties of the spectra of 
SNe 2001gh and 2002hy are significantly different from those 
of any other SN considered in Figure [3 

In Figure [6] we compare the I/-band absolute light 
curve of the same sample of objects shown in Figure [3 We 
may immediately note the large heterogeneity in these light 
curves. In particular, SN 200Igh is brighter than any other 
but SN 2002bj, and shows the slowest post peak decline. 


Based on the comparison described above, we consider 
SN 2001gh an unprecedented event. The object shares spec¬ 
troscopic similarity with SN 2002hy but their light curves 
are very different. If the spectra of the two SNe may sug¬ 
gest similar gas properties at the photosphere, their pho¬ 
tometric evolutions are also indicative that the progenitor 
parameters and/or the explosion energy may have been sig¬ 
nificantly different. In this context, we note the remarkable 
photometric similarity of SN 2002hy with the light curve 
of the Type Ibn SN LSQ13ccw (Figure (3), whose lumi¬ 
nous peak and fast photometric evolution were explained 
in ter ms of interaction betwe en the SN ejecta and He-rich 
CSM llPastorello et al.l I2OI5I I. Although SN 2002hy does 
not show the narrow He I emission from the CSM lines 
observed in Type Ibn SNe such as LSQ13ccw, we can¬ 
not rule out that a similar mechanism may also power the 
light curve of SN 2002hy. In this case, the lack of narrow 
lines from the photoionised ga s may be explained with a 
very high-density CSM ( e.g. see Blinnikov fc SorokinaluOlO : 
IChevalier Sz Irwinll201ll : iMoriva et al.ll2013f ). 

Is this scenario suitable also for SN 2001gh? SN 2001gh 
has a luminous and broad light curve peak which could in¬ 
dicate a moderate mass of ®®Ni, but the non-detection at 
~165 days after maximum light suggests that the ejected 
®®Ni mass cannot be the primary source powering its light 
curve. In fact, we will see in section [3 that the ejecta-CSM 
scenario suitable for Type Ibn SNe can provide a promising 


Figure 5. Comparison of the early optical spectrum of SN 2001gh 
with those of the peculiar Type Ib SNe 2002bj, 2002hy, 2005bf, 
and the Ibn LSQlSccw at epochs between -5 and 11 days from the 
maximum light. Vertical solid lines mark the rest frame position of 
the main He I features, while the dotted, dashed, and dot-dashed 
lines give the rest wavelength positions of the main H and Fe II 
features, along with those of the Call H&K doublet. All spectra 
have been reddening and redshift corrected assuming the values 
adopted in this work and from the literature. 



Figure 6. V-band absolute light curve for SN 2001gh compared 
with those of the peculiar SNe Ib 2002bj, 2002hy (H-band ab¬ 
solute magnitudes), 2005bf, and the Type Ibn LSQlSccw (sloan 
r-band converted into Johnson-Cousins i?-band). The magnitudes 
of SNe 2001gh and 2002hy have been corrected for the reddening 
values adopted in this paper, E(B — V) of 0.016 and 0.141 mag, 
respectively, and accounting for their adopted distance moduli of 
39.33 ± 0.10 and 33.61 ± 0.20 mag. Distances and extinction esti¬ 
mates for all other SNe have been taken from the literature. Ages 
are relative to the rest-frame B maximum light. 


explanation also for SN 2001gli-like events, although alter¬ 
native mechanisms will be discussed. 
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5 CONSTRAINING THE POWERING 

MECHANISM OF SN 2001GH 

It is clear from the data illustrated in the previous sections 
that SN 2001gh is a peculiar SN Ib with unprecedented ob¬ 
servational properties, showing some spectroscopic similar¬ 
ity with the usual SN Ib 2002hy. However, the latter has a 
much faster evolving light curve, similar to those of some 
Type Ibn SNe, s uch a s SN 1999cq, 200 0er and LSQlSccw 
(iMatheson et al.l l2000l : iPastorello et al.l l2008l . l2015ll , whilst 
SN 2001gh has a broad light curve wit h a slow post-peak de - 
cline, reminiscent of some SLSNe fe.g. iQuimbv et al]l201lll . 
From the discussion in sectional we suggested that the ob¬ 
served properties of SNe 2001gh and 2002hy are not primar¬ 
ily driven by the radioactive decays. A way to verify this 
claim is through data modelling. 

In order to constrain the parameters of the explosion 
of SN 2001gh and, hence, the progenitor’s scenario, we 
have estimated an approximated bolometric light curve 
for SN 2001gh. Since only two photometric bands were 
available for the light curve of SN 2001gh, the estimate of 
the bolometric luminosity required some assumptions and a 
few operational steps. Firstly, we derived the BV light curve 
by calculating the luminosity from the extinction-corrected 
apparent magnitudes in these bands and for our adopted 
distance. Then, since the SN 2001gh spectra are practically 
featureless, we computed the total luminosity from 900 to 
26000 A of two black bodies with the temperatures derived 
for each of the two spectra of SN 2001gh (see section 
[3dl . We estimated the bolometric corrections at these two 
epochs, and applied them to the BV psewdo-bolometric 
luminosities following a similar criteria than the one applied 
for the K-correction, i.e., we assumed the same correction 
derived from the first spectrum for the early epochs, the 
correction derived from the second spectrum for the late 
epochs, and for the intermediate epochs, we integrated the 
total flux of black bodies with intermediate temperatures 
assuming a smooth evolution of the continuum temperature. 
The errors in the bolometric luminosity account for the 
uncertainties in the estimates of the distance, extinction, 
apparent magnitudes and black body temperature. The 
inferred bolometric light curve is shown in Figure [T] The 
luminosity at maximum light was estimated by htting 
the light curve with low order polynomials, resulting in 
Tmax « (1.7±0.9)xl0‘‘® erg s-A 

Since the limited data set available does not allow us 
to robustly constrain the explosion scenario, we attempt to 
compute the explosion parameters and the nature of the 
progenitor star by considering three alternative scenarios for 
explaining the observables of SN 2001gh, in particular its 
bolometric light curve: 

• Radioactively-powered SN scenario — In our first 
model, we assume that the luminosity of SN 2001gh is mostly 
powered by the energy released in the radioactive decay 
chain ®®Ni - ®®Co - ®®Fe. In this case, we compare the light 
curve of SN 2001gh with the output light curves o btained 
with a toy model based on a simplified trea tment of lArnettI 
l|l982h and described in I Valenti et al.l l|2008l’F^ . This model, 

This model uses the analytic approximations of lArne t3 lll982ll . 


which has been successfully used for modelling a number o f 
stripped envelope SNe fe.g. iMorales-Garoffolo et al.ll2014l l. 
divides the SN evolution into a photospheric and nebular 
phase, assuming expanding ejecta with spherical symmetry. 
During the photospheric phase, the model includes the en¬ 
ergy produced by the radioactive decays. During the nebular 
phase, instead it assumes that the SN luminosity is powered 
only by the energy deposition of 7 -rays from the ®®Co de¬ 
cay. Accounting the photospheric velocities derived before 
(see section ITTI) . and the host galaxy distance and extinc¬ 
tion values as discussed in section we obtain a fairly 
good match with the bolometric light curvj^ adopting a 
kinetic energy between 0.6x10®^ and 3.5x10®^ erg, ejecta 
mass of about 2 Mq, and ®®Ni mass ~ 0.6 Mq (see panel 
a of Figure [TJ. Although this model reproduces fairly well 
the obtained bolometric light curve, the small ejecta mass, 
composed by about 30% in mass of nickel, does not easily 
accommodate in a typical CC-SN explosion scenario. Sim¬ 
ilar M(®®Ni)/Mejecta ratio can be found, in fact, in som e 
rare luminous Type la SNe le.g. iTaubenberger et al.ll2013l h 
However, in contrast with what expected in a thermonu¬ 
clear SN explosion scenario, there is no evidence of strong 
features of iron-group elements in the spectra of SN 2001gh. 
In addition, the total ejected mass here derived is compara- 
ble with the avera ge value for a normal Type Ib SN (see e.g. 
iDrout et al.ll20llh . However, to reproduce the broad lumi¬ 
nosity peak and the slow spectral evolution of SN 2001gh, 
one would have expected a larger ejected mass. In summary, 
for different reasons, the model illustrated above does not 
support a scenario where the SN 2001gh luminosity is pri¬ 
marily powered by radioactive decays. 

• Magnetar spin-down scenario — An alternative 
possibility to explain the properties of SN 2001gh is that the 
energy powering its luminosity comes from the spin-down of 
a newly formed magn etar, instead of the cano nical radioac¬ 
tive decays (see e.g. iKasen fc Bildst^ l201Clh . We explore 
this possibili ty using a m a gneta r mo del as presented an d 
discussed bv llnserra et al.l (l2013l l and iNicholl et aP ll2014h . 
The observed light curve of SN 2001gh is compared to a grid 
of magnetar-powered synthetic SN light curves, through 
minimization. The observed light curve satisfactorily fits a 
model obtained with ejected mass of ~ 3.4 Mq, a magnetic 
field B « 11x10^'* G, and a spinning period of Pi « 12 ms 
(see panel b of Figure [T]). This B field is high but in combi¬ 
nation with the Pi estimated, this magnetar could produc e 
a luminosity peak of 10"*® erg s“^ dKasen fc BildstCTl201Cll h 
which is in agreement with the maximum luminosity derived 
with our bolometric light curve. However, this scenario alone 
does not clearly explain the slow spectroscopic evolution of 
SN 2001gh or the boxy profiles of the He I lines. 

• Ejecta-CSM interaction scenario — The 
radioactively-powered models predict a faster evolu¬ 
tion of Vph, which is not observed in the available spectra of 
SN 2001gh. If we consider another observational constraint 
from the spectroscopy of SN 2001gh, i.e. some evidence of 
boxy profile in the Hel lines, an alternative interpretation 

acco unting for the typo correction reported in lArne^ dlQQdT) (see 
also IWheeler, Johnson &: Clocchiattill2015^ 

In order to fit the late phases of the SN, we considered in all 
cases the last upper detection limit as a genuine detection. This 
is intended to test the most conservative scenario. 
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could be offered, according to which SN 2001gh is a more 
canonical stripped envelope SN exploded in a dense, He-rich 
circumstellar shell. This CSM was produced through mass- 
loss episodes occurred prior to the SN explosion. To support 
this suggestion, we modelled the bolometric light curve of 
SN 2001gh using a semi-analytic code which accounts for 
ejecta-CSM collision as the main powering source (panel 


!Cta-UbM collision as tne mam powering source (panel 
of Figure [3 see Chatzopoulos. Wheeler fc Vinkol I 2 OI 2 I 


or iNicholl et al.l l2014l . and references thereinl. To obtain 
a good match between the light curve output of this 
interaction model and the observed peak of luminosity of 
SN 2001gh would require CSM/ejecta mass ratios in the 
range 0.1-0.2 to 1, assuming always a modest kinetic energy 
and an opacity of 0.2 cm^ g“^. For example, in Figure 
[7] we plotted two of the best interaction-powered models 
with ~ 7 Mq of ejecta and ~1 Mq of CSM (kinetic energy 
~ 0.25x10®^ erg), and with Mejecta = Mcsm « 2.4 M© 
(kinetic energy ~ 0.1x10®^ erg)- 


Therefore, accounting for the comparison between the 
observables of SN 2001gh (in particular its bolometric light 
curve), and the expectations from the models, we consider 
more plausible that the luminosity of this SN has been 
mostly powered by the interaction of the ejecta with an 
opaque CSM, probably ejected by the progenitor star years 
prior to the SN explosion. Accordingly, this scenario is sup¬ 
ported by the Hel line profiles observed in the SN 2001gh 
spectra. The low P-Cygni velocities and the modest velocity 
gradient estimated from the spectra would indicate that the 
lines formed in shocked CSM rather than in the SN ejecta. 
In summary, we cannot provide conclusive arguments to rule 
out some the alternative scenarios, in particular the magne- 
tar spin-down powering mechanism. On the other hand, the 
radioactive decays alone cannot provide a reasonable expla¬ 
nation for the properties of this peculiar Ib SNe. 


6 SUMMARY 

In this paper we have presented spectroscopic and photomet¬ 
ric observations of the peculiar SN Ib 2001gh, discovered at a 
redshift z=0.16 by the STRESS survey. Assuming a smooth 
evolution, SN 2001gh shows an unusual broad and delayed 
light curve peak, reaching —18.55 mag at B maximum light. 
On the other hand, the two spectra available for SN 2001gh, 
taken at around 30 days from each other, show a quasi fea¬ 
tureless blue continuum, with a few weak P-Cygni lines that 
we attribute to Hel at a velocity of 6000-7000 kms“^. 

Three possible explosion scenarios have been investi¬ 
gated to explain the observables of SN 2001gh: according 
to the first one, the SN luminosity may arise from radiated 
energy from the ®®Ni to ®®Co to ®®Fe decay chain, where the 
explosion ejected a total mass of ~ 2 M© (about 0.6 M© of 
which being ®®Ni); in the second scenario, the SN luminosity 
is powered by a magnetar with a magnetic field B ~ 11x10^^ 
G and a spinning period of about 12 ms; the third scenario 
suggests a more canonical SN Ib embedded in a dense, cir¬ 
cumstellar shell. Weighting up all the observables and the 
resulting models of SN 2001gh, and accounting for the spec¬ 
troscopic similarity with the peculiar SN 2002hy (and more 
marginally with Type Ibn SNe), we believe that the scenario 
with the peak SN luminosity being powered mostly by the 
interaction of the ejecta with an optically thick shell is the 


most reliable. Unfortunately, this dense CSM does not allow 
us to identify the true chemical composition of SN 2001gh, 
which is key information to unveil the nature of the progen¬ 
itor star. 

SN 2001gh is an unprecedented object. Assuming that 
this event is representative of a distinct class, we can derive 
an indicative estima te the rate of occurren ce by using the 
recipes described in iBotticella et ^ (l2008l b The main un¬ 
certainty in the calculation is the incomplete coverage of the 
light curve that we address by exploring two extreme cases 
for the evolution after maximum, eit her a decline simila r 
to SNe Il-Plateau such as SN 2004et ( Maguire_etal. gOld l. 
or to SNe II-Linear like SN 1979C (iBarbon et al.lll982l l. In 
both cases, the light curve is scaled to the observed magni¬ 
tude at maximum of SN 2001gh. With this prescription, and 
considering the observing log and detection efficiency of the 
STRESS survey, we derive a rate for SN 2001gh-like events 
to be in the range 1-4x10“® yr“^Mpc“®, which is 30 to 100 
times smaller that the overall rate of core collapse SNe at 
reshift z ~ 0.2. 

The outstanding properties of SN 2001gh are an indi¬ 
cation that our current knowledge on the mechanisms pro¬ 
ducing the variety of observed parameters of CC SNe is still 
incomplete. Object like SN 2001gh are intrinsically rare, and 
panoramic surveys and spectroscopic classification programs 
should help to discover further members of this family. 
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APPENDIX A: SN 2002HY 

SN 2002hy was discovered by iMonardI (120021 ') on 2002 
November 12.10 UT in NGC 3464. Its coordinates are: 
a = lO^'Sd^dOMS,)! = -21°03'41'.'2 (J2000.0; Figure [AT]|. 
The object was initially spectroscopically classified as a 
peculiar Type Ib a few days after the explosion. The 
spectrum, in fact, was dominated by a blue continuum 
with super-imposed stro ng He l lines (iBenetti et al.l l2002l : 
iHarutvunvan et al.ll2008h . We note that a few months prior 
to the discovery of SN 200 2hy, on 200 2 January 21.4 UT, 
another object (SN 2002J , IWeisz fc Lill200jl was observed 
in the same hos t galaxy, and cla ssified as a Type Ic near 
maximum light (ICappellardl2002h . 

NGC 3464 is a rather luminous (Ms ~ —21.4 mag) 
barred spiral galaxy, with a distance modulus of ^ = 33.61 ± 
0.20 ma£3- We also assume a total foreground extinction 
E{B — V) = 0.141 ± 0.030 mag, including a contribution 
form the Galactic reddening of E{B — V)g = 0.048 mag 
llSchlaflv fc Finkbeinerj [ 2 OI 1 I L and an additional contribu¬ 
tion of E{B — V)ngc 3464 = 0.093 inside the host galaxy. 
The latter was estimated from the equivalent width of the 
narrow Nam interstellar featur e and the relation given by 
IPoznanski. Prochaska fc BloomI (|2012| L 

To estimate the explosion date of SN 2002hy, we rely 
on the similarity of its J7-band light curv e with that of 
LSQ13ccw fFigure [6l IPastorello et Sll2015^ . Being the ex¬ 
plosion date of LSQlSccw well constrained, and assuming 
a similar rising time for the two SNe, we estimate the ex¬ 
plosion of SN 2002hy to occur on JD = 2452582.7 Ig (i-e., 

NED, NASA/IPAC Extragalactic Database; 

http://nedwww.ipac.caltech.edu/ 


SN 2002hy | 



Figure Al. /-band image of SN 2002hy field taken with LCO 40 
inch -I- CCD on 2002 November 18. The SN position is marked. 

2002 November 04.20 UT). This date is consistent with the 
discovery date on 2002 November 12 (JD = 2452590.6), and 
the last non-detection of the SN on 2002 October 13 (JD 
= 2452560.6). Deep images of SN 2002J taken on 2002 Jan¬ 
uary 25 with the ESO New Technology Telescope (NTT) and 
EMMI (ESO Multi-Mode Instrument) were also checked, 
showing no source detected at the SN position. 

One spectrum of SN 2002hy was taken at the ESO Very 
Large Telescope 3.6-m telescope with EFOSC2 (Faint Ob¬ 
ject Spectrograph and Camera - v.2; using grism 11 and 12) 
on 2002 November 15.33 (JD = 2452593.83; Figure [5ll, at 
phase ~ 11 days after the presumed explosion epoch. The 
spectrum was reduced following the prescriptions illustrated 
in section [ 33 ] 

The photometric data of SN 2002hy were reduced as 
detailed in section For the amateur data, an image of 
NGC 3464 taken on 2004 May 10 with a 0.3m telescope 
(-|- SBIG ST-7 CCD) was used as a template to subtract 
the host galaxy contamination. The apparent magnitudes 
of the SN in the R-band and the associated errors are 
reported in Table lAll According to the comparison in 
Figure |5] SN 2002hy seems a photometric twin SN of the 
Type Ibn LSQlSccw, having a fast average decline after 
maximum of 'ya 10.6 ± 0.5 mag/lOOd, comparable with 
tho se of the Type Ibn (b etween 8 and 13 mag/lOOd; e.g. 
see IPastorello et aI.ll201a L Therefore, considering also the 
information on the photometric evolution, we propose 
to re-classify SN 2002hy as a transitional Type Ibn/Ib 
SN, although with He I spectral lines which are broader 
than observed in other Type Ibn SNe (see e.g. Figure [5]). 
This sub-class of transitional events will be discussed in a 
forthcoming paper (Pastorello et al. in preparation). 
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Table Al. i?-band photometry of SN 2002hy. 


UT Date 

JD 

2400000.0+ 

Phase^ 

(days) 

R 

Instr. 

13/10/02 

52560.6 

-22.1 

< 18.0 

L.A.G.Monard 

12/11/02 

52590.6 

7.9 

16.28 + 0.18 

L.A.G.Monard 

14/11/02 

52592.6 

9.8 

16.66 ± 0.17 

L.A.G.Monard 

15/11/02/ 

52593.9 

11.1 

16.91 ± 0.13 

ES03.6m 

17/11/02 

52595.6 

12.8 

17.00 + 0.14 

L.A.G.Monard 

24/11/02 

52603.3 

20.5 

17.82 + 0.23 

AAT 

02/12/02 

52610.6 

27.8 

18.72 + 0.34 

L.A.G.Monard 

13/12/02 

52621.8 

39.0 

19.63 + 0.29 

J.B.Vanssay 


^ Relative to the explosion date (JD = 2452582.7). 
L.A.G.Monard = 0.3m Telescope + SBIG ST-7, 0.91^^pix“^; 
ES03.6m = ESO 3.6m Telescope -h EFOSC2, 0.33^^pix~^.; AAT 
= Anglo-Australian Telescope + WEI, 0.25"pix~^; J.B.Vanssay 
= 0.28m Telescope + KAF-400, 0.85^^pix“^ 
f Date in which the spectrum was taken. 
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